Construction of optimal witness for unknown two-qubit entanglement 
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Whether entanglement in a state can be detected, distilled, and quantified without full state recon- 
struction is a fundamental open problem. We demonstrate a new scheme encompassing these three 
tasks for arbitrary two-qubit entanglement, by constructing the optimal entanglement witness for 
polarization-entangled mixed-state photon pairs without full state reconstruction. With better effi- 
ciency than quantum state tomography, the entanglement is maximally distilled by newly developed 
tunable polarization filters, and quantified by the expectation value of the witness, which equals the 
concurrence. This scheme is extendible to multiqubit Greenberger-Horne-Zeilinger entanglement. 

PACS numbers: 03.65.Ud, 03.67.Mn, 42.65.Lm 



A fundamental issue is how to detect, distill, and quan- 
tify entanglement in a state [1], quantum correlation 
not imitated by classical correlation. These tasks are 
essential in understanding quantum nonlocality versus 
local realism [2], and in developing quantum technolo- 
gies where entanglement needs to be enhanced and mon- 
itored. An important open question is whether the tasks 
can be done without referring to the full knowledge of 
the state [3] . It is motivated from the demand of efficient 
protocols that directly access relevant quantities. 

Current schemes for those tasks rely on the full knowl- 
edge of the state estimated by quantum state tomography 
(QST) [4], even for two-qubit entanglement. For exam- 
ple, the detection by an entanglement witness [5], a Her- 
mitian operator whose expectation value is positive for 
all separable states [6.-^], requires QST to choose proper 
measurement settings. Procrustean distillation [ol, llO|, 
in which entanglement is enhanced by local filtering, has 
been also performed with known entanglement or QST. 
Since QST may have redundant information and becomes 
impractical for systems with larger number of qubits [sj , 
those approaches are indirect and inefficient. Hence, the 
schemes not relying on QST are necessary; there is a 
quantification scheme [ll', Ti] , which is however not al- 
ways available as it requires two state copies. 

It has been found [13l4l6| that entanglement detection 
and quantification can be done using a physical observ- 
able known as optimal witness Wp. Wp is defined relative 
to a state p as Tr(Wpp) = min^^^M Tr(Wp), where M is 
a collection of witnesses W and the trace Tr(-) gives the 
expectation value. For a two-qubit state p, it give s a pop- 
ular entanglement measure, concurrence C{p) |l4l Jl7Hl9| , 



C{p) = max{0, -2Tt{Wpp)}. 



(1) 



The proper collection M here contains all the witnesses 
which are obtained from the Bell witness WBeii (= — 
I Bell) (Bell I) through stochastic local operations and clas- 
sical communications (SLOCC), a composition of local 
unitary operation and local filtering. In Eq. ([1]), C{p) is 
interpreted as the "distance" between p and the support- 



FIG. 1. (Color Online) (a) Entanglement quantification by 
optimal witness Wp. Each witness W provides a supporting 
hyperplane Tt{Wp) = (solid lines) of the set of separa- 
ble states in the space (filled area) of two-qubit states. The 
"distance" d = max{0, —2Ti(Wp)} (thick arrows) between a 
state p and a hyperplane is maximized by Wp. The maximum 
distance equals C{p). (b) Construction of Wp via Procrustean 
distillation. The transformation from Bell witness Wseii to 
Wp is equivalent to that from p po) to the maximally dis- 
tilled state pdis (= P2), the state optimally detected by Wseii. 
In our scheme, the distillation is achieved iteratively using 
tunable local filters; the subscript is the iteration index. 



ing hyperplane TiiWpp') = of the set of separable states 
in the space of two-qubit states; see Fig. [TJa). One can 
obtain C{p) by constructing Wp, rather than indirectly 
from a calculation based on QST estimation. 

It has not been known how to construct Wp without 
QST. We will do it for two-qubit entanglement, based on 
the following link between Wp and Procrustean distilla- 
tion |^,[iQ|; see Fig.djb). The optimal witness Wp^.^ rel- 
ative to pdis, the maximally distilled state obtained from 
p, is constructed from WBeii by local unitary operation. 
Moreover TiiWpp) and Tr(Wp^.^pdis) are identical, up to 
constant 5o; see below. These facts, originated from the 
SLOCC invariance of C, mean that the distillation oper- 
ation connecting p and pdis is equivalent to the filtering 
part of Wp. We obtain pdis to construct Wp below. 

In this Letter, we provide a new entanglement detec- 
tion scheme not relying on QST, based on Wp and indi- 
vidual copies of a state. For a unknown two-qubit mixed 
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FIG. 2. (Color Online) Setup. In the source, polarization- 
entangled photon pairs are generated from spontaneous para- 
metric down- conversion (SPDC) in two cascaded beta-barium 
borate (BBO) crystals pumped by pre-walk-off-compensated 
pulses. After spatial-mode filtering by single-mode fibers 
(SMFs), decoherers convert them into a mixed state p. Then, 
entanglement in p is distilled by local filters, in which a tun- 
able polarization-dependent loss can be applied with an arbi- 
trary basis using beam-displacing prisms (BDPs). Finally, 
filter-transmitted photon pairs are analyzed with quarter- 
wave plates (QWPs), half- wave plates (HWPs), polarizing 
beam splitters (PBSs), and SMF-coupled photon counters. 
Our scheme never uses the QST estimation of p (inset). 



state p of polarization-entangled photon pairs, we con- 
struct Wp and VFpdis (hence determine C{p) and C{pd[s))^ 
by maximally distilling p and then measuring two-qubit 
correlation with applying local unitary operation. This 
approach is more efficient than QST for the distillation 
and quantification, and extendible to multiqubits. 

We describe the experiment setup in Fig.O The source 
generates an unknown mixed state p of polarization- 
entangled photon pairs, in which qubit |0) (|1)) is en- 
coded in horizontal H (vertical V) photon polarization. 
Here, a pure Bell state of photon pairs, prepared through 
SPDC [23, [in, becomes p due to asymmetry of pump 
polarization and severe decoherence applied in decoher- 
ers [21] composed of birefringent quartz crystals (22j . 
Next, photon j (= 1,2) of p moves to local filter j. 

Our scheme requires that any arbitrary filtering 
(SLOCC) operation is realized by the local filters. For 
this purpose, we have newly designed the filters [23]. By 
filter jf, photons having polarization parallel to its filter- 



ing basis are filtered off (not proceeding to detector j) 
with probability 1 — Pj, while those with orthogonal po- 
larization are always transmitted. The filtering basis and 
transmission probability pj of filter j are tunable over the 
whole range by a QWP and a HWP and by two BDPs 
and the HWP in between, respectively, while preserving 
interferometric stability; see the inset of Fig. [3l^a). 

In the detectors, filter-transmitted photon pairs un- 
dergo local unitary rotation of photon polarization, re- 
alized with a QWP and a HWP, and are projected onto 
a two-qubit product state by PBSs. Coincidence count 
(^1^2) is obtained from photon number nj detected at 
counter j. For our severely decohered state p, we choose 
the measurement time of 10 s (during which the source 
creates 5 x 10^ photon pairs) per detector setting. 

We first maximally distill p iteratively, based on the 
fact [24] that the degree of polarization (DOP) of each 
photon of pdis vanishes, i.e., all the local reduced density 
matrices are fully mixed. We first (at step k = 1) erase 
DOP 1 of photon 1 by tuning filter 1, next (at k = 2) do 
the same for photon 2, and then (at /c = 3, • • • ) repeat- 
edly re-adjust filters 1 and 2 alternately, until both the 
DOPs vanish within experimental uncertainty [25^]. At 
each step for photon j, the filtering basis and transmis- 
sion probability pj for erasing are determined from the 
measurement of the local reduced density matrix [22]. 
This iterative distillation finishes usually within a few 
steps and does not require the full reconstruction of p. 

After the distillation, Wp is constructed with local uni- 
tary operation on photon polarization. We consider the 
two-qubit correlation of (4(nin2) — 2(ni) — 2(712) +A^)/A^, 
which has the information of two-photon Stokes param- 
eters (ni) ((^2)) is the number of photons detected 
at detector 1 (2) and N is the total number of photon 
pairs passing through the filters among M pairs injected 
from the source. There appear three different pairs of lo- 
cal extremum values of the two-qubit correlation, ±Ai,2,3 
(Ai > A2 > A3 > 0), if one measures it with varying lo- 
cal unitary operation (by tuning the detector QWPs and 
HWPs). At the step of the maximal distillation. 



So 



Tr(I^pP) = ^(l-Ai-A2+gA3), 



(2) 



where sq = N/ {M^pip2) and the sign factor q = —I 
for all entangled states [25]. Moreover, Tr(Wp^.^pdis) = 
(1 — Ai — A2 -\-qXs)/4: = Tt{Wpp)/so. Hence, by measur- 
ing A's, one constructs both Wp and Wp^.^ and determines 
C{p) and C(pd\s) = C{p)/so using Eq. ([T]). We measure A's 
efficiently [25| , by first determining the detector settings 
for A's and then measuring A's. We perform 16 inde- 
pendent detector settings of coincidence counting for the 
former, and 12 settings (4 settings/A x three A's) for the 
latter. The 12 settings for A's are determined from the 
16-setting data through the singular value decomposition 
of two-photon Stokes parameters [22] . Crude determina- 
tion of the settings for A's is enough for precise detection 
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of A's, since A's are extremum values; a small error (~ 6) 
of the setting determination causes a much smaller error 
(~ (^^) of the detection of A's. 

Before discussing our experimental data, we intro- 
duce the concept of the best witness. Often one can 
explore not the whole set M of witnesses but only 
its subset due to certain limitation, and may 
not include the optimal witness Wp. In this case, 
we propose to use the best witness Wp,M' within M', 
Tr(VFp,M'P) = minvFGM'Tr(VKp). For two qubits, it gives, 
via max{0, — 2Tr(Wp^M'P)}, the best lower bound of C(p), 
which is the best among those achievable by the wit- 
nesses in M^ In our experiment, the witness constructed 
via Eq. (j2j) is the best witness at the steps where the 
maximal distillation is not achieved [22|, and it becomes 
the optimal witness at the step of the maximal distilla- 
tion; see below. This supports that our determination 
of C{p) does not exceed the exact value, even if there is 
experimental imperfection. And, it allows us to extend 
our scheme to the quantification of three-qubit entangle- 
ment, for which it has not been known how to achieve 
Wp. Note that the new theoretical findings of this work 
include the link between Wp and the scheme in Ref. 
the best witness, and the extension to multiqubits [22 

We discuss our experimental results. DOPs 1 and 2 
approach zero for larger k; see Fig.[3l They vanish below 
< 0.1 already at k = 3,4, implying that the maximal 
distillation is almost done. This conclusion is supported 
by the fact [25] that the error caused by partial distilla- 
tion with DOPs < 0.1 is of the order of (DOP)^ in the 
determination of C{p) by Wk- It is also supported by 
the pattern of (nin2) around the local extrema A's [see 
Fig. 13]: At /c = 0, classical correlation of local qubit in- 
formation leads to a complicated pattern, while at /c = 3, 
it is removed, resulting in a simple pattern. 

At each step /c, we measure A's and 5o, and plot 
TT{Wkp) = (1 - Ai - A2 + ^A3)/5o and TT{Wdis,kPdis,k) = 
TT{Wkp)/so in Fig.lHb). We show that -2TT{Wkp) pro- 
vides a lower bound of C{p) [22|. The lower bound is the 
largest among those obtained from the witnesses con- 
structed with varying local unitary operation under the 
filtering setting fixed at step k. Namely, Wk and Wdis,/c 
are the best witnesses at step k. At the maximal dis- 
tillation step, Wk = Wp. For comparison, Cs are es- 
timated [18] by QST and maximum likelihood estima- 
tion at each k. The distillation enhances concurrence 
from C{p) 0.18 to C(pdis) ^ 0.33, and -2TT{Wkp) 
approaches to the QST-concurrence C{p) within 0.03 at 
k = 3,4. Given experimental uncertainty (see below), 
the optimal witnesses Wp and Wp^.^ are constructed at 
/c = 3,4, and directly provide C{p) and C(pdis)- 

We mention the imperfections in the data, (i) the in- 
complete erasure (< 0.1) of DOPs, (ii) the difference 
(< 0.03) between the concurrence by Wk=3^4 and that by 
QST, and (iii) the fiuctuation (< 0.03) of'c(p) by QST 
at different /c's. They might come from unwanted group 
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FIG. 3. (Color Online) Distillation and Quantification, (a) 
Degree of polarization (DOP) and filter transmission proba- 
bility pj (inset) at distillation step k. At /c = 0, no filtering 
is applied, (b) The values of -2TT{Wkp), -2Tr(Vl/dis,fcPdis,fc), 
C(p), and C(pdis,fc); Wk and Wdis,fe are the best witnesses for 
p and the distilled state pdis,k at step k, respectively, and Cs 
are estimated by QST. Error bars represent statistical uncer- 
tainty [di (counts) ""^^^l. Data points in (b) has the same size 
as the error bar (±0.004) for VK's at /c = 0. For larger k, the 
error bars for VK's increase due to smaller N/M. 



delay [23| (< 10 /im, depending on the filter parameters) 
between the polarization components in the local filters 
and/or random fiuctuation of coincidence counts due to 
the Poissonian statistics; other error sources may exist. 
Imperfection (i) might enhance (ii), as it leads small error 
of ~ (DOP)^ in the determination of C{p) by Wp] this 
error is not included in the error bars in Fig.[3jb). Given 
these uncertainties, it is reasonable to stop the distilla- 
tion at the step {k = 3) of DOPs < 0.1. Our Monte Carlo 
simulation reproduces qualitatively the same features. 

We discuss the efficiency of our scheme. At k = 3,4, 
the error bar of C{pdis) determined by Wp^.^ is about three 
times smaller than that (- ± 0.02) by QST. This is ob- 
tained with the time cost of 520 seconds in our scheme 
(240 for three distillation steps, 160 for determining the 
settings for A's, 120 for measuring A's) and 160 seconds 
in QST (only for quantification of pdis)- Our scheme is 
more efficient than QST, as its error bar will be roughly 
two times smaller than that of QST with 520 seconds. 
Our simulation indeed shows [22] that for the distilla- 
tion and quantification of pdis, QST requires a number of 
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We have shown how to efficiently identify unknown 
two-qubit entanglement. It is the ffist time, to our knowl- 
edge, that the optimal witness has been constructed for 
entanglement quantification without QST. 

We thank J. Eisert, P. Kwiat, L. C. Kwek, Tzu-Chieh 
Wei, and H. Weinfurter for discussions, NRF (Mid-career 
Researcher, 2009-2332), and KRISS project ("Single- 
Quantum-Based Metrology in Nanoscale"). H.S.P. and 
S.S.B.L. equally contributed to this work. 



FIG. 4. Two-photon coincidence pattern. Coincidence counts 
of (nin2) in a subspace (^i, O2) of the parameters of local uni- 
tary operation. The subspace is chosen to contain the posi- 
tions of two of the local extremum values of the two-qubit 
correlation such that A^, — A^, A^/, and — A^/ are located at 
{61,02) = (0,0), (7r,7r), (7r/2,7r/2), and (37r/2, 37r/2), respec- 
tively: (a),(d) / = 1 and = 2, (b),(e) / = 3 and = 1, (c),(f) 
/ = 2 and = 3. Three panels (a)-(c) are for step /c = 0, while 
(d)-(f) for /c = 3. At step /c = 3 with maximal distillation, 
(nin2) has the same information as the two-qubit correlation. 
Panels (d)-(f) show a simple symmetric pattern, clearly dis- 
tinct from (a)-(c), and reveal the quantum correlation in p 
masked by classical correlation of local qubit information. 

measurements roughly two times larger than our scheme. 
The better efficiency comes from the fact that in our 
scheme the relevant quantities (A's, which linearly de- 
pend on measured counts) are directly accessed, using 
Wp and the SLOCC invariance of C, while in QST the 
estimation of C requires redundant information such as 
local properties of p and is sensitive to the estimation 
error due to the nonlinear dependence of C on p. 

Our scheme is extendible to Greenberger-Horne- 
Zeilinger (GHZ) entanglement [26|; see Ref. [22]. The 
iterative distillation works for A/'-qubit (pure or mixed) 
states [24. It is useful for detection of 'hidden' non- 
locality [9[ of multiqubits. The quantification is ex- 
tendible to three-qubit GHZ entanglement; it does not 
work for W-state entanglement, which behaves differently 
from GHZ under SLOCC. For a pure state, the opti- 
mal witness is obtained from a GHZ witness VFqhz = 
3//4 - I GHZ) (GHZ I by applying the SLOCC operation 
for maximal distillation, then the exact value of a GHZ 
entanglement measure Ts ^] is obtained without QST. 
For mixed states, the same procedure gives a reasonably 
good lower bound of T3. This quantification is powerful, 
as it is a formidable task to compute a GHZ entanglement 
measure for mixed states even with the QST information. 
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